Bifurcations in Boltzmann-Langevin One Body dynamics for fermionic systems 
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Wc investigate the occurrence of bifurcations in the dynamical trajectories depicting nuclear 
collisions at Fermi energies. The quantitative description of the reaction dynamics is obtained 
within a novel one-body transport approach, based on the solution of the Boltzmann-Langevin 
equation in three dimensions, with a broad applicability for dissipative fermionic dynamics. 

Dilute systems formed in central collisions are shown to fluctuate between two energetically 
favourable mechanisms: reverting to a compact shape or rather disintegrating into several frag- 
ments. The latter result can be connected to the recent observation of bimodal distributions for 
quantities characterising fragmentation processes. 



Phase transitions are general phenomena occurring in 
interacting many-body systems [1-4]. Over the past 
years, many efforts have boon devoted to the identifi- 
cation of new features related to finite-size effects. As 
shown by recent thermodynamical analyses, first-order 
phase transitions in finite systems are characterised by 
negative specific heat and bimodal behaviour of the dis- 
tribution of the order parameter [5, 6]. The latter physi- 
cally corresponds to the simultaneous presence of differ- 
ent classes of physical states for the same value of the 
system conditions that trigger the transition (like the 
temperature, for instance). 

In particular, the appearance of phase transitions from 
the liquid to the vapour phases has been widely investi- 
gated in the context of the nuclear multifragmentation 
phenomenon [4, 7-9]. Indeed, due to the analogies be- 
tween the nuclear forces and the Van-der-Waals interac- 
tion, the nuclear matter equation of state (EOS) fore- 
sees such a possibility [10, 11]. The theoretical find- 
ings cited above have stimulated corresponding thermo- 
dynamical analyses of the properties of the products is- 
sued from nuclear reactions at Fermi energies. Under 
suitable conditions, a bimodal character of experimen- 
tal observables, such as the size of the heaviest cluster 
produced in each collision event [12], or the asymmetry 
between the charges of the two heaviest reaction prod- 
ucts [13] has been revealed. Many investigations have 
also been focused on the complex nuclear many-body 
dynamics, to probe the reaction mechanisms governing 
the occurrence of phase transitions [14 17]. Within such 
a context, nuclear fragmentation studies have recently 
pointed out that bimodality could have a dynamical ori- 
gin, related to the; fragment formation mechanism [18], 
without necessarily requiring the reaching of thermody- 
namical equilibrium. 

From a general point of view, interacting many-body 
systems may experience a very rich dynamics, ranging 
from mechanisms dominated by one-body (mean-field) 
effects to phenomena governed by strong fluctuations and 
correlations. In the regime of low-energy collective pro- 
cesses, nuclear dynamics presents a rather stable charac- 
ter; this is the domain where the fluctuation mechanism 



can be described in the small amplitude limit, restrict- 
ing to mean-fleld (quantum) fluctuations of collective ob- 
servables [19, 20]. This limit is exceeded when violent 
perturbations, like for instance dissipative heavy-ion col- 
lisions, bring the system beyond the one-body collective 
dynamics, with two-body nucleon collisions and correla- 
tions playing an important role. Along the compression- 
expansion path traced by the nuclear reaction, fluctua- 
tions introduce the anisotropy seeds from which 'nuclear 
droplets' can develop. More precisely, the system may 
access mechanically unstable regions of the EOS, called 
spinodal, where a density rise is related to a pressure fall; 
there, phase-space fluctuations are even amplified, lead- 
ing to phase separation [8, 17]. As soon as a mottling 
pattern stands out at low density, i.e. at the boundary 
of the phase separation, a bundle of bifurcations into a 
variety of different dynamical paths may set in. 

The aim of this work is to further investigate the 
dynamical trajectory of disassembling nuclear systems, 
seeking for features associated with phase transitions. 
Thus we will explore the possible occurrence of bifur- 
cation patterns and bimodal behaviour in central heavy- 
ion reactions at Fermi energies. This study is undertaken 
in the framework of a novel transport approach to deal 
with the dynamics of quantum many-body systems, the 
Boltzmann-Langevin-One-Body (BLOB) model. Start- 
ing from the procedure introduced in [21], the BLOB 
model represents a numerical solution in three dimen- 
sions of the full Boltzmann-Langevin (BL) equation. The 
latter is a transport equation that describes the time evo- 
lution of the semiclassical one-body distribution function 
/(r, p, t) in response to the mean-field potential, includ- 
ing fluctuations and correlations due to hard two-body 
scattering [22-25]. 

Hence the distribution function / evolves according to 
the action of the effective Hamiltonian ff[/], the aver- 
age Boltzmann collision integral /[/], and the fluctuating 
term 8I[f] as: 

f = d,f-{H[f]J} = I[f] + 6I[f]. (1) 

Like in standard transport approaches [26] , we sample 
the dynamics through the test-particle method, under the 
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assumption of spatial and temporal locality of the two- 
body coUisional process. Attest particles per nucleon are 
employed. As firstly proposed by Bauer and Bertsch [27], 
in order to solve the BL equation, nucleon-nucleon (N- 
N) correlations are included by organising test particles 
in phase-space agglomerates of iVtcst elements of iden- 
tical isospin which simulate nucleon wave packets; this 
arrangement is redefined at successive time steps and lo- 
cally, for couples of colliding agglomerates. According 
to this rescaling and for elastic N-N collisions only, the 
average rate of change of the occupancy fa around the 
phase-space location {ra,Pa) at a given time takes the 
form: 



/a(rQ,Pa) = 9 



dpb 



dil W^(AB«-i.cD) _F(ab^cd) 



(2) 

where g denotes the degeneracy, and integrations are over 
momenta pb and scattering angles fl. The first integra- 
tion argument is the symmetric transition rate from an 
intermediate state AB to a final configuration CD. To se- 
lect the test particles defining the nucleon wave packets 
A and B, we adopt the following procedure: We consider 
a sphere, centred at the position r^, with radius equal to 
the scattering distance, associated with the free elastic 
N-N cross section at Fermi energies (taken equal to 50 
mb); among all test particles inside the sphere, we pick 
up the A'tost closest particles to the elements a and b in 
momentum space, respectively. The final state is rep- 
resented by c e C and d G D. The transition rate is 
obtained by averaging over all couples of test particles 
involved in the transition T, = (AB — ^ CD), in terms of 
relative velocity and differential N-N cross section: 



all test particles belonging to the agglomerates A and B 
can be reconsidered as starting points of new collision 
processes, the scattering probability has to be suitably 
rescaled, dividing it by N^^^^. Once the sorting allows for 
a scattering to occur, modulation functions are applied 
to precisely adapt the density profile of final-states to 
the available vacancy profile /a...d, with the requirement 
of imposing the most compact configuration compatible 
with the constraint of energy conservation [29]; the re- 
sulting occupation functions of the modulated final-state 
density profiles (/a...d)m should approach the unit. The 
accurate treatment of Pauli-blocking factors represents 
the major improvement of the present approach with re- 
spect to the method proposed in [27]. We notice that 
this is an essential point for keeping a correct description 
of the dynamics of fermionic systems [30] , as also shown 
in [21]. 

The extension of the wave packets makes necessary to 
pay special attention to scatterings close to the surface 
of the system, i.e. occurring across potential boundaries: 
confronting the shape of the wave packet to the shape of 
the surface, the blocking factors are increased in propor- 
tion to the spread of the nucleon packet outside of the 
boundary (similarly to what is done in some molecular- 
dynamics approaches [14]). 

In the following, the BLOB model is applied to a well 
documented but highly constraining phenomenology: the 
low-energy threshold for multifragmentation in head-on 
heavy-ion collisions at Fermi energies. We will also refer 
to results of the so-called Stochastic Mean Field (SMF) 
model [17, 25], that corresponds to an approximate treat- 
ment of the BL equation. In SMF, fluctuations arc pro- 



VF(AB<-i-CD) — 



.da 



\ = 



(3) 



The second integration argument contains the product of 
occupancies of the entire agglomerates /a...d and of the 
associated vacancies /a...d: 



F(aB^Cd) = /a/b/c/d-/a/b/c/d = (F(ab^cd) 



(4) 

Rewritten in terms of test-particles, the representa- 
tion of eq. 4 indicates that only the fraction of the pack- 
ets which are really modified by the scattering can sig- 
nificantly contribute to the transition probability, while 
overlapping volumes contribute to the Pauli-blocking fac- 
tors. 

On this basis, full phase-space fluctuations are intro- 
duced in the equation of motion by moving simultane- 
ously the test-particle agglomerates, in analogy with the 
extended TDHF procedure of including perturbations in 
the Slater conflguration [28]. The scattering is decided 
by confronting the probability W x F with a random 
number and scanning the entire phase space in search of 
collision configurations at successive time steps. Since 
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FIG. 1: (Color online) N-N collisions in "''Xe-H^^-'Sn at 
32^MeV for b — 0: mean Pauli-blocking rejection rate (Rp), 
and mean evolution of the number of effective collisions (Nch) 
which equals its variance cr^ (A^ofi). Bottom band. Time evo- 
lution of the projected density profile. Right band. For the 
same initial conditions, bifurcations lead to several different 
fragmentation patterns at 300 fm/c. 
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jected on the coordinate space and injected by agitating 
the spatial density profile. We will show that the im- 
plementation of fluctuations in full phase space (as in 
BLOB) improves significantly the reaction dynamics. 

Like in ref. [31], for the simulation we use a soft equa- 
tion of state with fcinf = 200 MeV; the potential compo- 
nent of the symmetry energy of the EOS is given by a 
linear term as a function of the density (asy-stifi'). For 
numeric purposes, the transition rate has been rewritten 
as ~ {\va — Vb\)s{das/dfl) by introducing a screened 
in- medium cross section as (from ref. [32] ) , corresponding 
to the transition (AB o CD) . 

Fig. 1 illustrates a study of the N-N collision statistics 
for the collision i36xe-)-i24gjj ^t 32 AMeV, at a central 
impact parameter (6 = 0). The number of events consid- 
ered is 600. Only a tiny fraction of the attempted N-N 
collisions is permitted by Pauli-blocking factors. Indeed 
the Pauli-blocking rejection rate Rp, defined as the ra- 
tio between rejected and attempted collisions, is close to 
the unit. The average number of effective N-N collisions 
(TVeff) reflects the average density evolution of the sys- 
tem: a maximum is encountered in correspondence with 
the largest compression and a minimum marks the subse- 
quent expansion mechanism. This evolution is illustrated 
for one single event in the lower band of fig. 1. Later 
on, (Nas) increases again indicating that the process of 
fragment formation sets in, and it finally levels off till, 
around 300 fm/c, the system has stabilised into a well 
defined configuration. Thus fragments are formed when 
the system expands, according to the spinodal decom- 
position scenario [17]. The equality between the mean 
value and the variance ensures that Langevin fluctuations 
stand out with the correct amplitude [21]. Some exam- 
ples of fragmentation patterns appearing at 300 fm/c for 
the same initial conditions are collected in the right band 
of the flgure; their variety already signs the presence of 
bifurcations. 

Having verified that Langevin fluctuations show up 
with the expected properties, we let vary the incident 
energy over a large interval, from 19 to 58AMeV. The 
corresponding evolution of the fragment multiplicity, cal- 
culated by counting all fragments with element number 
Z > 4: standing out at 300 fm/c, is shown in flg. 2. The 
general behaviour that BLOB traces is that a transition 
from incomplete fusion (just one nucleus observed) to 
multifragmentation occurs between 20 and 30AMeV in 
i36Xe-h^^^Sn. The less neutron rich system '^^'^Xe+'^'^'^Sn, 
also simulated, presents the same evolution, but retarded 
with a shift of about 2 ^MeV to larger bombarding en- 
ergies. Indeed the corresponding IMF multiplicity is 
smaller, in agreement with experimental findings [33]. 
The same observable is also deduced at asymptotic time, 
by letting the system cool down after 300 fm/c (where the 
excitation is around 3.5MeV per nucleon) through a pro- 
cess of sequential evaporation, simulated in-fiight in the 
swarm of light ejectiles and fragments (for this purpose. 
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FIG. 2: (Color online) Top. Multiplicity distribution of 
fragments with Z > 4 as a function of the incident en- 
ergy in ""^Xe-F^^^Sn {b = 0) at 300 fm/c (filled contours) 
and after secondary decay (contour lines), Bottom. Corre- 
sponding mean values; additional calculation for the system 
^■^*Xe-|-^^^Sn is also shown. The heaviest system is also sim- 
ulated within the SMF approach. 

the model SIMON [34] was used). In the cold system 
the fragment multiplicity results increased by secondary 
decay up to around 45AMeV and levels off for larger 
bombarding energies due to a more prominent decay into 
light fragments. The system evolves from a dominance 
of incomplete fusion around 20^MeV to a full multi- 
fragmentation pattern. This result is in close accordance 
with experimental observations addressed to the same 
system, where multifragmentation is already seen above 
30^MeV [33]. If the BLOB approach is replaced by a 
SMF calculation with corresponding mean-field param- 
eters, the transition moves above 40AMeV. Hence the 
formation of inhomogeneities is enhanced in the BLOB 
approach and, beside a possible additional amplification 
from spinodal modes, this process is induced by fluctua- 
tions which are implemented in full phase space. More- 
over, fluctuations have the effect of anticipating the for- 
mation of fragments. Indeed, the fraction of energy spent 
in light-particle emission reduces in favour of a larger con- 
tribution to the developing of inhomogeneities and of a 
more explosive dynamics, which drives the separation of 
those inhomogeneities into fragments. On the contrary, 
for small bombarding energies (below 30AMeV), due to 
insufficient kinetic energy, the mean field often succeeds 
to coalesce the inhomogeneities into a compact system. 
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Fluctuations feed the large width of the multiplicity 
distribution above the multifragmentation threshold. A 
better insight on their distribution, and on the possi- 
ble mixing of different reaction mechanisms is provided 
by the charge asymmetry a = (Zi — 2-2)1 [Z^ + .^2) be- 
tween the largest fragment Z\ and the second largest 
fragment Zi emitted in each event. This observable, pre- 
sented in fig. 3 illustrates that the reaction mechanism 
changes from incomplete fusion (charge asymmetry close 
to the unit) to fragmentation (very small charge asym- 
metry) discontinuously, passing through a bombarding- 
energy range (around 28 AMeV ) where both mechanisms 
coexist. Thus we observe, at this transition energy, a 
wide, but discontinuous, bunch of trajectories and a bi- 
modal behaviour of a, indicating that either the system 
recompacts or it breaks up into several pieces of similar 
size. Intermediate charge- asymmetry values are never 
populated. It should be noted that this variety of con- 
figurations is related to the same "macroscopic" initial 
conditions, i.e. the same beam energy and impact pa- 
rameter. These features, often discussed in the context 
of the thermodynamics of phase transitions and bimodal 
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FIG. 3: (Color online) Top. Survey of the reaction mecha- 
nism: distribution of charge asymmetry a as a function of 
the incident energy in "''Xe-l-^^^Sn (6 = 0) at 300 fm/c. Bot- 
tom. Corresponding mean values. Additional calculation for 
the system '^^'*Xe-|-^^^Sn, and simulation of the heaviest sys- 
tem within the SMF approach and after secondary decay are 
also displayed. 
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FIG. 4: (Color online) Temporal evolution of the size distri- 
bution of potential concavities in ^^^Xe-|-^^''Sn (6 = 0) for an 
incident energy of 32 ^MeV (filled contours) and 25 AMeV 
(contour lines). 

behaviour of the order parameter [6, 12, 13], are observed 
here as a result of the fragmentation dynamics, governed 
by the spinodal decomposition mechanism. As a con- 
sequence of the presence of fluctuations, the amount of 
light-particles early emitted may vary from one event to 
the other, leading to energy fluctuations in the fragment- 
ing system. Then the latter behaves as in a thermal bath 
and, at the transition energy, oscillates between two con- 
figurations which are energetically favourable. The same 
behaviour is also found in the corresponding SMF ap- 
proach, but at larger incident energy. After secondary 
decay, the transition pattern is not washed out if events 
producing large fission residues are removed, and it could 
be searched in experimental data relative to the same sys- 
tems [33]. 

To explore in finer details the developing of instabili- 
ties, the evolution of the size of concavities of the mean- 
field potential is followed in time for the incident energies 
32 and 25 AMeV. The distribution of the corresponding 
radius is shown in fig. 4. This radius probes the size 
of forming blobs of matter which are the fragment nest- 
ing sites; they may eventually separate into clusters and 
leave the system (prominent mechanism at 32AMeV), 
or merge into larger fragments, or fuse back together 
forming a large residue (as clearly seen at 25AMeV). 
The potential concavity-radius evolves initially with the 
composite-system size, as marked by the upper branch in 
fig. 4, and suddenly drops to smaller sizes which charac- 
terise all forming droplets of matter within a small vari- 
ance, as expected for a spinodal multifragmentation pic- 
ture: this latter defines the lowest branch in fig. 4, and 
reflects in the bimodal pattern of the charge-asymmetry 
observable shown in flg. 3. 

In conclusion, a promising new framework of transport 
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modelling has been realised by extending the Boltzmann- 
Vlasov formalism to include a Langcvin residual term 
acting in full phase space and adapted to fermionic sys- 
tems. This approach reveals to keep the specific char- 
acter of a one-body theory, as far as the description of 
mean-field (spinodal) instabilities in concerned, and at 
the same time to solve the long-standing difficulty of in- 
troducing fluctuations of correct amplitude in a Boltz- 
mann formalisation. The model, tested on the transition 
from fusion to fragmentation in central collisions at Fermi 
energies, reveals to be closer to the observation than pre- 
vious attempts to include a Langevin term in Boltzmann 
theories. Moreover, we have identified the occurrence of 
bifurcations and bimodal behaviour in dynamical trajec- 
tories, linked to the fragment formation mechanism. At 
the transition energy, the system may either recompact 
or split into several pieces of similar sizes. We would like 
to mention that our findings do not challenge the validity 
of thermodynamical analyses evidencing the occurrence 
of bimodal behaviour in nuclear fragmentation [12, 13]. 
Indeed a large fraction of the available phase space is 
populated through the spinodal decomposition mecha- 
nism [8], thus legitimating the use of thermodynamical 
equilibrium concepts. Our calculations provide a possi- 
ble explanation, based on the occurrence of dynamical 
instabilities, of the origin of trajectory bifurcation and 
bimodality. Finally, we stress that the results presented 
here may be relevant not only for nuclear fragmentation 
studies, but in general for the dynamical description of 
quantum many-body systems. 

Clarifying discussions with M.-F. Rivet, B. Borderie 
and J. Rizzo are gratefully acknowledged. 
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